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It is shown that for collision cascades the global fractal dimension cannot give an adequate de- 
scription of the geometrical structure because it is insensitive to the internal anisotropy of the object 
arising from the directionality of cascade branches. In order to give a more elaborate description 
of the cascade we introduce an angular correlation function, which takes into account the direction 
of the local growth of the branches of the cascades. It is demonstrated that the angular correlation 
function gives a quantitative description of the directionality and the interrelation of branches of 
the cascade. The power law decay of the angular correlation is evidenced and characterized by an 
exponent j3 and an angular correlation length R a different from the radius of gyration R. It is 
demonstrated that the overlapping of subcascades has a strong effect on the angular correlation. 

PACS numbers: 47.53. +n, 47.54. +r, 61.43.+Hv 



I. MOTIVATIONS 

Collision cascades develop in condensed matter as a 
consequence of irradiation with energetic beams of par- 
ticles. The bombarding particles transfer their kinetic 
energy in series of collisions with the target atoms and 
the energized, recoiling atoms generate further recoils in 
their own slowing-down process. The result of this energy 
sharing process is a collision cascade. The cascade can 
also be considered as a far-from equilibrium process lead- 
ing to structure formation in the solid target. Recently 
the geometrical structure of collision cascades has been 
analyzed by means of analytical approaches and Monte 
Carlo simulations concentrating on the possible fractal 
and multifractal aspects. These investigations have been 
extended to the study of the self-similarity properties 
of the cascade to the determination of its fractal 

dimension and multifractal spectra for different interac- 
tion potentials and they gave an insight into the 
cascade-subcascade transition and into the spike creation 

Ml- 

Far from-equilibrium processes often create complex 
geometrical structures which exhibit fractal properties 
characterized by the fractal dimension D JlQ] . However, 
the fractal dimension is a global property of a cluster of 
particles, it does not provide an insight into the struc- 
tural details of the object. It has been demonstrated in 
the case of Diffusion Limited Aggregation (DLA) that 
large branching structures may be internally anisotropic. 
This anisotropy results in a tangential correlation dif- 
ferent from the radial one fllf| and it shows up in the 
behavior of the three point correlation function as well 
& 

The internal structure of collision cascades, the cor- 



relation of the vacancies, the interrelation of branches 
and the directionality of ion tracks can have an impor- 
tant impact on several physical processes, e.g. ion-beam 
mixing |l3) , radiation enhanced diffusion jl3],[l4) , electric 
and optical properties of irradiated polymeric materials 
p5| , ^6[ , the development of mechanical stress due to ion 
bombardment and crack propagation in the solid from 
the irradiated zone [ p7{ . 

In the case of diffusion, for example, the correlation 
of vacancies due to the directionality of the cascade 
branches can enhance the diffusion process in the irra- 
diated zone. Inside the cascade the path of the displaced 
atoms can be viewed as conduit, forming a labyrinth 
in the damaged material jL4|. Thus, we are led to de 
Gennes's model of "termite diffusion" , when the diffusing 
particles (termites) meet the conduit-labyrinth, they dif- 
fuse much more quickly and over much greater distances 
than they do in the bulk JTsJ . 

Further motivations to study the directionality of ion 
tracks in a solid comes from the earlier studies of system 
of particles interacting via elastic collisions. A system 
where many particles scatter, interacting elastically with 
each other, may be regarded as a model of gas. In such a 
system each particle is scattered successively, so it walks 
almost randomly. Although the velocity autocorrelation 
of such a particle was believed to decay exponentially, 
Alder and Wainwright discovered a long time tail by com- 
puter simulation Jl9| , p0[ : 

(v(t) ■ v(0)) ~ r d ' 2 , d>2, (1) 

where d is the dimension of the embedding space. Long 
time tails of random walks are also found in a much sim- 
pler system called the Lorentz gas model where a classi- 
cal particle is scattered elastically off randomly located 
fix scatterers. (This model was originally introduced as a 
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model of electron motion in a metal.) In the case where 
the scatterers form a regular lattice, it is proved that 
the velocity autocorrelation has exponential decay. On 
the contrary if the scatterers are distributed at random 
forming a fractal with dimension D there is a long time 
tail for the negative values of the scalar product of the 
velocities Q: 

(v(t) -v{0)) ~ (2) 

Hence the memory effect in the directionality seems to 
be a common feature of the motion of particles scattered 
randomly. 

Based on the above motivations, in the present paper 
we want to obtain a more elaborate picture of collision 
cascades gaining information about the internal branch- 
ing morphology. The computer simulations have the ad- 
vantage that one can monitor such quantities which are 
hard to measure or are not measurable but they can help 
the understanding of a physical process. In the study of 
ion-solid interaction in certain cases not only the distri- 
bution of the vacancies is of interest but even the path 
of the moving particles. Obviously, the study of the 
ions's paths includes the vacancies as well, the position 
of the vacancies are the branching points in the cascade 
tree. Thus in the present study not only the vacancies 
but all the collisional points are included in the analysis 
of the cascade regardless of the vacancy creation, that 
can easily be done by means of computer simulations. 
We performed Monte Carlo simulations of self ion col- 
lision cascades in two and three dimensions. For the 
three-dimensional simulations the scattering cross sec- 
tion da = C(m)E~ m T~ 1 ~ m dT belonging to the inverse 
power law potential V(r) ~ r -' v l m was applied. Here E 
is the energy of the bombarding particle, T is the energy 
transferred to the target atom in the scattering process 
and m denotes the parameter of the interaction potential 
< to < 1. In two dimensions we established a toy model 
with the scattering cross section da — KE^~T^ 1_ ~ dT 
0. A two dimensional cascade does not have experi- 
mental relevance but it can help to reveal the role of the 
embedding dimension d in the formation of the cascade 
structure. For further details about computer simulation 
of ion-solid interaction see also Ref . 0| . 

In the following, at first the structure of cascades gen- 
erated by Monte Carlo simulations is analyzed in terms 
of the density-density correlation function. Then an an- 
gular correlation function is introduced based on the di- 
rection of the local growth of the cascade branches. It 
is shown that the angular correlation function gives a 
quantitative description of the directionality and the in- 
terrelation of branches of cascades, furthermore, the role 
of the intersection of branches is enlightened. 



II. CORRELATION FUNCTIONS 

The usual way of studying the structure of an object 
composed of N particles in the d dimensional embedding 
space is in terms of the density- density correlation func- 
tion C(r): 

°M = W^gt E P(ri)p(rj), (3) 

where p(r) is 1 if there is particle at r and zero otherwise. 
In our case these particles are the points where collisions 
occurred in the target material during the evolution of 
the cascade, fid denotes the solid angle in d dimension. 
For structures of finite size, like collision cascades, the 
two point density-density correlation function depends 
on the overall size of the structure described by a charac- 
teristic length R as well as the internal length r. In this 
case the correlation function can be written as C (r, R) . 
Here C(r, R) represents the correlation function averaged 
over a large number of cascades of the same size R. The 
characteristic macroscopic length of the cascade is the 
radius of gyration R defined as the average distance of 
the particles: 

* 2 = 4£kW/~^/^C(r) (4) 

id 

C(r, R) is normalized by the total number of particles N: 

N = fl d J drr d - x C{r) (5) 

Assuming that C (r, R) is a homogeneous function of its 
variables results in the scaling form: 

C(r,R)=r- a g(x), x = ~, (6) 

where g(x) is constant for x « 1 and g(x) ~ e~ x for 
x > 1. The function g{x) is called scaling function and 
the exponent a is the scaling exponent. The fractal di- 
mension D of the object is defined through the behavior 
of the total number of particles N as a function of the 
radius of gyration R, i.e. N ~ R D . From Eqs. (||||) it 
follows that a = d — D. If Eq. (^) provides a valid de- 
scription of the geometric scaling properties then plots of 
r a C(r, R) vs. r/R for structures of different sizes will fall 
on a common curve (on the scaling function g(x)). This 
data collapse provides a reliable way of measuring the 
scaling exponent a (and thus the fractal dimension D) 
by varying the value of a until the best collapse is ob- 
tained. Representative examples of the density-density 
correlation function and its data collapse analysis can be 
seen in Figs. |l|, [2] for the 2d and 3d models. The ex- 
cellent collapse obtained with ten different system sizes 
proves the validity of the scaling ansatz of Eq. (^) for 
both model systems. 
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However, this way of description based on C(r, R) 
cannot reveal anything about the possible internal 
anisotropy because the density-density correlation func- 
tion is insensitive to the structural details. In order 
to obtain information about the branching structure we 
have to take into account some dynamical features of the 
growth of the cascade. For this purpose we assign to each 
particle, i.e. to each point of collision, with position fi 
the unit vector p(f) of the linear momentum of the scat- 
tered particle. If there is vacancy creation in a collision 
the unit vector of the recoiled particle appears at the po- 
sition of its first collision. With this prescription one and 
only one unit vector is assigned to each point of the ob- 
ject studied. The p(f) vector characterizes the direction 
of the local growth of the corresponding branch of the 
cascade tree at position f during the cascade evolution. 
Fig. |H shows the unit vectors p attached to the parti- 
cles. Using the vector field p(fi) an angular correlation 
function C a (r, R) can be introduced with the following 
definition: 



Ca(r) 



1 



Nrd-ifldSr 



E 

-<\ri-rj\<r+± 



p(n)-p{f 3 ) (7) 



The form of this equation is similar to Eq. (j^) but in the 
summation the scalar product of the p(f) vectors is used 
instead of the product of the one particle densities p(fi). 

Along a given branch of the cascade the vectors p(f) 
are correlated in the sense that they are almost parallel 
to each other hence the scalar product in ([?]) has values 
close to 1. Subbranches appear in the cascade tree as 
a result of the recoiled particles. In Fig. || it can be 
observed that a high energy recoil can give rise even to 
an extended subcascade the branches of which grow in- 
dependently of the other parts of the cascade, the other 
subcascades. This implies that the average value of the 
scalar product of vectors p{fi) belonging to different sub- 
cascades is close to zero or it can even be negative. In this 
sense the angular correlation function C a (r, R) can mea- 
sure the directionality of cascade branches with respect 
to each other. Representative examples of the absolute 
value of the angular correlation function \C a (r, R) \ of cas- 
cades generated in the 2d model are shown in Fig. [|. Fig. 
|]a presents a comparison of the angular and the density- 
density correlation functions belonging to the same sys- 
tem and Fig. [p> shows \C a (r, R)\ for different system 
sizes R. It can be observed that for a certain length range 
C a {r, R) decays according to a power law with an expo- 
nent significantly different from the exponent describing 
the algebraic decay of C (r, R) . At a certain value of r 
the angular correlation function becomes negative and it 
goes to zero through negative values at long distances. 
Let's denote R a the characteristic length where C a (r,R) 
becomes negative. Increasing the system size R in Fig |]b 
R a is also increasing. The negative part of C a (r, R) has 
a local minimum (i.e. \C a {r, R)\ has a local maximum) 



and it goes to zero exponentially due to the finite size 
cutoff. 

Based on the above observations, the same form of scal- 
ing ansatz can be assumed for C a (r, R) as for C(r, R), i.e. 
Eq. (§): 



C a {r,R) = r-^f{x), 



r 

R' 



(8) 



where f{x) denotes the scaling function and a a the scal- 
ing exponent belonging to C a (r, R) . At small distances 
C a (r,R) approaches C(r, R) thus part of the power law 
decay of C a (r, R) can be attributed to the power law de- 
cay of C(r,R), therefore the exponent (3 characterizing 
the angular correlation is the difference of a a and a: 



C a (r,R) 
C(r,R) 



(3 = a a 



a 



(9) 



Finally, the internal anisotropy of cascades is determined 
by two quantities, by the characteristic length R a where 
C a {r, R) drops to negative values and the exponent (3 
describing the speed of the decay of C a (r, R) . 

In the next section we present the results of the anal- 
ysis of systematic Monte Carlo simulations in terms of 
the correlation functions introduced above, and we try 
to give a simple interpretation of the results obtained. 



III. RESULTS AND DISCUSSION 

In a recent publication (see Ref . ) it has been demon- 
strated that the parameter to of the interaction poten- 
tial plays the role of a control parameter from the view- 
point of the cascade geometry. For decreasing m a struc- 
tural transition takes place in the cascade from an open 
branching structure to a space filling one. When the up- 
per critical dimension d u of the cascade exceeds the di- 
mension of the embedding space d, so called geometrical 
correlations arise in the cascade due to the intersection 
of different branches and the overlapping of subcascades. 
These geometrical correlations lead to a fractal dimen- 
sion D(m) different from the self similarity dimension 
D (m) for to < 1 and to < i in the 2d and 3d models, 
respectively 0. The question naturally arises how the 
intersection of branches influences the behavior of the 
angular correlation function introduced. 

The above analysis in terms of the density-density and 
angular correlation functions was performed for cascades 
in a wide range of the parameter to of the potential ( 
in the same range as in Ref. Q]) for both model systems. 
The scaling exponents a, a a and the characteristic length 
R a were obtained by means of the data collapse method. 
Since the results concerning to the fractal dimension of 
cascades have already been discussed in Ref. M here we 
restrict ourself to the discussion of the properties of the 
angular correlation function. 
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Fig. |5| shows representative examples of the data col- 
lapse analysis of the angular correlation function in the 
2c? and 3d models. It can be observed that C (r, R) ex- 
hibits power law behavior for more than one order of 
magnitude in length. The good quality of the collapse 
obtained demonstrates the validity of the scaling ansatz 
Eq. (||). Note, that the positive and negative parts of 
C a (r, R) have the same scaling properties. The power law 
behavior of C a (r, R) expresses the fact that there is long 
range angular correlation inside the cascade, which man- 
ifests in the directionality of the branches. The negative 
regime of C a (r, R) can be interpreted as anti-correlation 
of the local growth directions at large distances. This im- 
plies that the angle of the unit vectors p(ri), p(rj) tends 
to be larger than ^ for \ fi — fj\ > R a . The characteristic 
length R a where C a (r, R) drops to negative values can 
be considered as an angular correlation length. Let us 
denote 7 the ratio of the angular correlation length R a 
to the radius of gyration R. 

From the scaling analysis, the exponent /3 characteriz- 
ing the decay of C (r, R) and the ratio 7 of the angular 
correlation length R a to the radius of gyration R can be 
determined precisely. From the simulations it was found 
that 7 is between 0.7 and 0.8 for both model systems 
and within the accuracy of the calculations we could not 
reveal any systematic dependence of 7 on the parameter 
of the potential to. 

Fig. |^ presents the values of (3 as a function of to for 
2d and 3d. It is important to note that small value of 
(3 corresponds to slow decay of the directionality, i.e. to 
long range angular correlation, while increasing j3 implies 
the weakening of the angular correlation. In Fig. ^ one 
can observe that in 3d two different regimes of (3 can be 
distinguished. If to > 1/3, when no overlap occurs in the 
cascade, (3 practically coincides with the self similarity 
dimension D a . In this region (3 is increasing for decreas- 
ing to. When m is close to one (Coulomb-scattering) 
the probability of creating a high energy recoil is rather 
small, therefore dense subcascades cannot occur, the cas- 
cade is strongly directed giving rise to a small value of (3. 
Decreasing m gives rise to an increasing number of sub- 
cascades growing independently, which results in weak- 
ening of the directionality order. This is expressed by 
the increasing value of (3. At to = 1/3 the self similarity 
dimension D Q and (3 reaches d/2, half of the dimension of 
the embedding space. For m < 1/3 when the overlap of 
the branches dominates the cascade structure, (3 remains 
practically independent of m, f3 = d/2. But this does 
not entails that the density-density correlation function 
C(r, R) and the angular correlation function C a (r, R) are 
independent of to. Their exponents a a and a are de- 
creasing with decreasing to but in such a way that their 
difference remains constant: 

= a a -a=^. (10) 



This argument is also supported by the two dimensional 
simulations (see Fig. ^a). It has been shown in Ref. 
that in the two dimensional model the overlap effect, the 
intersection of branches is always present, in the whole 
region of to. In Fig. [(^ it can be seen that (3 is equal to 
1, i.e. (3 — d/2, as it can be expected from the 3d case, 
see Fig. ||&. The result that (3 — | throughout the over- 
lapping regime demonstrates that the p{fi) vector field 
does not become isotropic even for high degree of overlap, 
i.e. for small values of to. In the system there is always 
a preferred direction, namely, the initial direction of the 
bombarding particle. This also manifests in the overall 
shape of the cascade. The overall shape is not spherical 
but it is always elongated along the initial direction of 
the bombarding particle. 

It has been mentioned in the introductory part that 
the internal anisotropy of DLA clusters was studied by 
means of the tangential and the three-point correlation 
functions, which are essentially special type of density- 
density correlation functions. It is important to mention 
that, beside these techniques, Meakin introduced an an- 
gular correlation function of the type of Eq. (0) p3[ . 
He assigned a so-called bond vector to each particle in 
the aggregate. The bond vector 6 is a unit vector point- 
ing in the direction of the bond, which is formed when a 
new particle is added to the growing cluster. Using the 
b vectors an angular correlation function could be de- 
fined. (3 — 0.31 was found in two dimensional off-lattice 
simulations |23| that evidenced the slow decay of the di- 
rectionality of the cluster branches in DLA processes. 

Collision cascades have the advantage that the struc- 
ture formation can be easily controlled by varying the 
parameter of the interaction potential to. This gives 
us a rich spectrum of possibilities to study the inter- 
nal anisotropy resulted from the directionality of the 
branches of the object and it makes also possible to gain 
information about the role of the intersection of branches. 

IV. CONCLUSIONS 

In the present paper we studied the internal anisotropy 
of collision cascades arising from the branching structure. 
It was demonstrated that the global fractal dimension 
cannot account for the internal details of the structure. 
To give a more elaborate description of the cascade it 
was necessary to introduce an angular correlation func- 
tion, which takes into account the direction of the lo- 
cal growth of the cascade branches. With the help of 
the angular correlation function we could quantify the 
anisotropy, which manifest in the directionality of the 
branches. The power law decay of the angular correla- 
tion was evidenced and characterized by an exponent (3 
and an angular correlation length R a different from the 
radius of gyration R. The overlapping of subcascades 
has a strong effect on the angular correlation. In the 
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absence of overlap (5 coincides with the self similarity di- 
mension, while in the presence of overlap (3 is constant, 
i.e. (3 — d/2. It is interesting to note that in the presence 
of overlap solely the dimension of the embedding space 
determines the speed of the decay of the directionality 
order. We argued that this internal anisotropy of growth 
directions shows up in the overall shape of the cascade 
as well. 
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VI. FIGURES 





FIG. 1. The density-density correlation function for system sizes from N — 400 to N — 4000 in the 2d model. The parameter 
m of the scattering cross section was chosen to be m = 0.7. The data collapse analysis was performed with ten different curves. 
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FIG. 2. The density-density correlation function for system sizes from N = 800 to N = 5000 in the 3d model. The 
parameter m of the scattering cross section was chosen to be m = 0.4. The data collapse analysis was performed with ten 
different curves. 




FIG. 3. Collision cascade generated in 2d. The arrows represent the unit vector p(r*;) of the linear momentum of the 
scattered particles after the collision, m has the same value as in Fig. [j]. The circle indicates an extended subcascade resulted 
from a high energy recoil. The surface of the solid and the direction of the penetration of the bombarding particle are also 
shown. 
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FIG. 4. Representative examples of the behaviour of the angular correlation function C a (r,R) in the same system as in 
Fig. [j]. (a) Comparison of the absolute value of the angular correlation function \C a (r,R)\ and the density-density correlation 
function C(r,R) belonging to the same system size R, (b) \C a (r,R)\ for different system sizes. 




FIG. 5. Data collapse analyzis of the angular correlation function in the 2d and 3d models. One can observe the good 
quality of the collapse, which demonstrates the validity of the scaling ansatz Eq. ^ . The parameter m of the scattering cross 
section was chosen to be m = 0.7 and m = 0.4 in 2d and 3d, respectively. The value of 7 = R a /R is indicated by dotted lines 
in the figures, it is 0.78 and 0.76 in 2d and 3d. 
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FIG. 6. The exponent /3 — a a — a characterizing the decay of the angular correlation function C a (r,R) (a) for 2d and (fe) 
for 3d. In the 3d case the curve of the self similarity dimension is also shown. The horizontal solid lines indicate the value of 
d/2. 
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